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,-H . Abstract. The status of small x resummation in the timelike kinematics is discussed. We present a 
^ ' general procedure to extract the large logarithms of x in the MS factorization scheme and to resum 
^^ ' them in a closed form. New results for the doubly-logarithm-resummed coefficient functions will 
^O , be reviewed. All our resummation formulae are in agrement with the fixed NNLO computations 
■ recently done by other groups in the MS scheme. 
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r-| ' To fix the ideas, we consider the cross section for the semi-inclusive hadron produc- 

'""'■ tion in electron-positron annihilation: 

"> '• e+ih) + e-{k2) ^ V\q) ^ h{ph) +X, (1) 

cni ; 

'sf . where V* is a virtual vector boson with virtuality Q = q = (^i + ^2) and X stands 

^ \ for any allowed hadronic final state. Here we are interested in the differential cross 

• ■ section for the single hadron production do^{x, Q^)/dx where x is the scaled momentum 

Q ■ fraction of the produced hadron h: 

Tl: x=^^, Q<x<\. (2) 

/\ • According to QCD facotrization the cross section can be written as the convolution of the 

j^ . partonic cross section to produce a parton k with scaled momentum fraction z = x/z' with 

the fragmentation function Dj^n^ {z!) from the parton the hadron h with scaled momentum 
fraction x (see e.g. Ref.[l]): 

Usually the partonic cross section ddg+^-^i^/dx is written in terms of the coefficient 
functions Ck{x, Q^) defined as 

c.(^'e) = -(^^^^, (4) 

where a^ contains all the electro- weak over-all factors. 



Perturbation theory fails when the fraction x of available energy carried away by 
the observed particle is too low, because large logarithms spoil the convergence of 
the perturbative series. The largest logarithms, the double logarithms (DLs), in the 
splitting functions that determine the evolution of the fragmentation functions have been 
computed to all orders a long time ago [2], and have even been used to perform LO 
global fits in QCD [3, 4] to data measured at the smallest x values. The DLs appearing 
in the coefficient functions Q(jc, Q^) are not resummed at LO and are expected to be not 
as important as those appearing in the evolution because they only appear at and beyond 
NLO in QCD. However, the inclusion of the DLs in the coefficient functions could 
make a significant improvement to the accuracy of cross section calculations making 
the analysis of Refs. [3, 4] feasible also at NLO. The complete DL contribution to 
partonic cross sections has been calculated in Ref.[5] for the case in which the coUinear 
singularities are regularized by giving a small mass nig to the gluon, the so-called 
massive gluon (MG) regularization scheme. The inconsistency noted in [1, 6] between 
the NNLO DLs calculated from the resummed result in Ref. [5] and those calculated 
from the fixed order result in Refs. [7, 8, 9, 10, 11, 12] is not surprising because the two 
computations were carried out in two different regularization and factorization schemes, 
namely the MG and the MS scheme. The DLs of the gluon coefficient functions have 
been computed for the first time in the MS scheme very recently by us in [6]. Here we 
will rederive the same result in a less formal but more simple way. 

Our goal is to get a resummed cross section in the MS scheme for which the dimen- 
sional regularization is necessary. To extract the leading logarithmic behavior, we exploit 
the factorization of the single gluon probability emission in the soft-coUinear limit. This 
is a consequence of the eikonal approximation and color coherence as it has been proven 
a long time ago in [13, 14]. This result can be used [6] to obtain the probability gluon 
emission in J = 4 — 2e dimensions: 

ll^\^ {AkY dx dz 
^) r(l-e);ci+2e? 



dw{x,z,e)^2%a4^^] '^.,^,, -^, (5) 



where as — Ofj/2;r, jl is the dimensional regularization scale taken here and in the 
following equal to the renormalization scale and where z = (1 — cos0)/2 with 6 the 
scattering angles of the emitted soft gluon with respect to the hard jet direction. Here 
% = Ca for a gluon jet and "^ = Cp for a quark jet. The expression for the probability 
emission given in Eq.(5) is what we need to obtain the gluon probability density in 
dimensional regularization. Fig.l shows a diagrammatic derivation of a consistency 
relation for the differencial cross section for gluon jet production which is: 

dOg = dOg +dOg^^dw{x,z,£). (6) 

Now introducing G{x,£) the gluon distribution density and taking the limit n — > oo, we 
obtain immediately the following bootstrap equation for it: 

x^+^'W{x,z,£) = 5(1 -x)+ f dx' f dz'K{x',z',£)x^+^'W{x',z\£), (7) 

Jx Jz 

where K{x, z, e) = dw{x, z, £)/dxdz and where the factor oix^^^^ represents our normal- 
ization coming from the explicit computation for the first gluon emission with n = 2. In 
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FIGURE 1. Here da" represents the cross section for the production of a gluon up to real corrections 
of order n, dPSt is the /-particles phase space and dw{PS\) ~ dw{x,z,£) after that azimutal integration is 
performed. The factorization of the single gluon probability emission has been used in the third step. 



Eq.(7) the value ^(x, z = 0, e) = G{x, e) should be taken only at the end of the computa- 
tion. This ensures the strong angular ordering of the emitted gluons, which is necessary 
to extract correctly the leading logarithms as proven in Refs.[2, 13]. Performing the 
Mellin transform, 

f((0)= I dxx'^fix); (0 = N-\, (8) 

of Eq.(7), then performing also the z-integrals, solving recursevely for ^(jc,z, e) and 
finally putting z = we obtain: 



G(a),e) = l + £ 



k=\ 



IUsCa 
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(9) 



According to the QCD factorization theorem we have that all the coUinear singularities 
in Eq.(9) should be factorized. In the MS factorization scheme this is done requiring 
Eq.(9) to be compared with [15, 16], 



G{co,£) = G^^ ( (0,a,, =j ) exp 



1 r^ii^li^FfSe da 



eJo 



f^ico, 



(10) 



where 5e = (4;r)^e^^^^ with Ye the Euler number and where jAf is the arbitrary factor- 
ization scale. The direct comparison of the two Eqs.(9,10) is non trivial. As shown in 
[6] a possible way to do this is to compare a simple differential equation satisfied by 



G(ft), e) given in Eq.(9) with the same differential equation this time obtained by use of 
the factoriation constraint given in Eq.(lO). We report here the result which is given by: 
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where 
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(12) 

Expanding the result in Eq.(ll) up to NNLO we obtain perfect agreement with the 
leading logarithmic terms of the fixed order result computed in the literature in the same 
scheme (see e.g. Eqs.(A.3,A.6) in Ref.[12]). 

We conclude noting that attaching the gluon jet to the quark lines of the LO process 
in Eq.(l), we have that the gluon coefficient function is according to Eq.(4): 
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(13) 



This result enables us to resum all the DLs in the gluon coefficient function in the MS 
factorization scheme for the first time and according to Refs.[17, 18] to have full control 
over all the large logarithms in all the coefficient functions at NLO. Our result Eq.(13) 
is also a key ingredient in the determination of the single logarithms in the timelike 
splitting functions, which is left to a forthcoming paper [19]. 
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